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Recently much progress has been made in applying field theory methods, first developed to study 
X-ray edge singularities, to interacting one dimensional systems in order to include band curva- 
ture elTects and study edge singularities at arbitrary momentum. Finding experimental confirma- 
tions of this theory remains an open challenge. Here we point out that spin chains with uniform 
Dzyaloshinskii-Moriya (DM) interactions provide an opportunity to test these theories since these 
interactions may be exactly eliminated by a gauge transformation which shifts the momentum. 
However, this requires an extension of these X-ray edge methods to the transverse spectral function 
of the XXZ spin chain in a magnetic field, which we provide. 



2 



I. INTRODUCTION 



One dimensional (ID) interacting systems exhibit unusual correlation effects dominated by strong quantum fluc- 
tuations. Fortunately, an array of powerful theoretical methods exist to study this physics, which is finding many 
experimental realizations. One powerful method is based on bosonizatioufi"— leading to the Luttinger liquid concept. 
Traditionally these methods are based on low energy effective field theory and only apply to the low energy excitations 
occurring near certain wave- vectors (such as q = 0). In the case of fermion models they begin by linearizing the dis- 
persion relation near the Fermi energy and ignoring irrelevant band curvature effects. However, in the last few years, 
these bosonization methods have been significantly extended by using techniques first developed to study X-ray edge 
singularities!^"— This has shown that band curvature effects, while formally irrelevant in the renormalization group 
sense,— can nonetheless have important effects on line-shapes of spectral functions even at low energy. Perhaps even 
more importantly, by combining these techniques with Bethe ansatz methods, it has become possible to make exact 
predictions of critical exponents at arbitrary momentumj^i^ That is, spectral functions are predicted to have the form 
S{q,uj) — >■ A[uj — a;i,(q)]~'' near singular energies wl(<z) for arbitrary q where both aji(q), which is not small, and the 
exponents r](q) are determined exactly using the Bethe ansatz. Recently, this approach has been extended to also 
obtain the amplitudes of correlation functions. A, using the Bethe ansatz»ii This new approach has been applied to a 
number of systems including fermions and bosons moving in the continuum ^^i^i^ii^ the fermion spectral function for a 
tight binding model^i and the longitudinal spectral function of the xxz S=l/2 spin chain in a magnetic field^ii^ with 
Hamiltonian: 

N 
i=l 

The longitudinal spectral function is: 

/oo 
^^^_iq,+i^t ^ s]{t)S^o{Q) > (1.2) 

3 

Wc set the lattice spacing a — 1. While this is accessible, at arbitrary q, to neutron scattering experiments these 
would normally measure a sum of S"^^ and the transverse spectral functions: 

/oo 
^^^~^qJ+^u^t ^ S]{t)S'„ (0) > (1.3) 

J 

for (s,s') = {+,—) and (—,+). For non-zero h, and S ^ are different. Due to the "Jordan- Wigner string 

operator" , the transverse spectral function is much more difficult to compute. While the high neutron fluxes available 
at the Spallation Neutron Source may eventually make such experimental confirmation possible, here we explore 
another route. 

Quasi ID magnetic compounds that do not have link parity symmetry (reflection about the midpoint of a link) will 
generally have magnetic Hamiltonians containing anti-symmetric Dzyaloshinskii-Moriy a^^'"'^^ interactions: 

6H = J2Dr{S,xS,+^). (1.4) 

i 

The two standard cases are staggered Dj = {—lyB and uniform Dj = D. Staggered DM interactions are invariant 
under site-parity (reflection about a site) but violate symmetry of translation by one site. On the other hand, uniform 
DM interactions have no parity symmetry whatsoever but respect full translation invariance. We restrict our further 
discussion to the case where Dj cx z so that the DM vector is parallel to the easy (or hard) axis of the symmetric 
exchange interactions, a situation which is sometimes dictated by symmetry. Then it is possible to eliminate the 
DM interactions exactly by a gauge transformation, yielding the standard xxz model of Eg. (jl.ip with modified 
parameters. This approach was used to study staggered DM interactions in Ref. ([l^ and [TtI ) and the theory was 
applied to a number of real materials. Here we consider the case of uniform DM interactions, Dj = zD and choose 
Z? > 0. In this case the gauge transformation is: 

5+=e-*"^5+ 5| = 5| (1.5) 

where 

a = tan"i(^). (1.6) 
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By use of above unitary transformations we have 

H^JY. (^^^^+^ + SfSf+, + ^eff~St~SU,) - h J2 St- (1.7) 

i i 

where the new exchange couphng and anisotropy parameters are given by 

J = ^/.P + D"^ Ae// = Acos(a). (1.8) 

The electron spin resonance (ESR) adsorption intensity, in standard Faraday configuration, is proportional to the 
transverse spectral function at q = 0, since the wave-vector of microwave photons is much less than the inverse 
lattice spacing. After the gauge transformation, the ESR intensity is therefore proportional to and S for 
the Hamiltonian of Eq. (|1.1|) at <? = a. [By using circularly polarized microwave radiation both 5''' (a, lo) and 
S ^{a,oj) could be measured separately] Thus the edge singularities predicted by X-ray edge methods at a non- 
zero wave-vector a given by Eq. (jl.6p are directly measured by ESR. ESR on spin chain compounds with uniform 
DM interactions therefore would provide a powerful probe of the new bosonization predictions. Quasi-ID spin- 
1/2 antifcrromagnctic insulators containing DM interactions with a uniform component include Cs2CuCl 4^^'2'^ and 
KCuGaFgjSi This provides a strong motivation to extend the X-ray edge methods to study edge singularities in the 
transverse spectral functions of the xxz chain in a magnetic field. 

In the next section we review results on the transverse spectral function using standard bosonization and then show 
that band curvature effects (in the equivalent fermion model) render these results invalid close to edge singularities. 
In Sec. Ill we apply X-ray edge methods to the model obtaining new results on the leading edge singularities. In Sec. 
IV sub-dominant singularities are discussed. Section V discusses ESR with uniform DM interactions, based partly on 
the results of Sec. III. Sec. VI contains conclusions and open questions. 



II. SPECTRAL FUNCTION OF THE XXZ SPIN CHAIN 



The xxz 5 = 1/2 model of Eq. (|l.ip is equivalent to an interacting spinless fermion model by the Jordan- Wigner 
transformation: 



Sf = c c 



1 



2 

57 = (-l)^exp(*7r^4cfc)c, (2.1) 
k<j 

The Hamiltonian (jl.ip is transformed to 

H = -iJ2i4c^+l + h.c.) -hY^c\c, 

+ JA^(clc,-l)(4iC,+i-i). (2.2) 

Note the factor of (—1)-' in the second line of Eq. (|2.ip . necessary for the first term in the fermionic Hamiltonian, Eq. 
(|2.2p to have the standard minus sign, h is the chemical potential of the fermionic model with h = corresponding 
to half-filling. For the non-interacting case A = the Hamiltonian is just a free fermion model and by going to 
momentum space 



, '1'- 



Co-e 



where p = lixnjN for periodic boundary conditions, the energy spectrum of non-interacting fermion model is found 
to be 

i? = ^(-Jcos(39)-/i)4cp (2.3) 
For this model the longitudinal spectral function S~''' is known exactly and the result is 



(2 J si 



sm S ~ LO'- 
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We see that S^^ is non-zero only in finite frequency range for a fixed momentum q, these lower and upper thresholds 
for kp < 7r/2 are given by 

uju{q) = 2 J sm — sm(fcF + — ) 

Wi(<7) = 2JsinMsin(fc^_M) (2.5) 

and for kp = n/2 

u)u{q) = 2 J sin -y 
uJhiq) = Jsin|g| (2.6) 

From Eq. ()2.4p . it is easily seen that the support of S"^^, for A = 0, is restricted to the interval [ujl^uju]. The nature 
of these lower and upper thresholds is different for zero and non-zero magnetic field, for non-zero magnetic field and 
kp < 7''/2 the lower threshold is given by creating one deep hole with momentum q, whereas the upper threshold 
is achieved by excitation of one high-energy electron with momentum q. But for zero magnetic field, kp — 7r/2, 
particle-hole symmetric case, the lower threshold is achieved either by creating one deep hole or one high energy 
electron but higher threshold is given by symmetric particle hole excitation around Fermi point with momentum of 
each q/2. 

In contrast to the longitudinal spectral function, even for A = and having the entire spectrum at hand, finding 
the transverse spectral function is a very difficult problem. The reason for this complication is the complicated form 
of the transverse spin operators, which have the string operator in their fermionic form Eq. (|2.ip . The equal time 
transverse correlation function, in this case, is known exactly22 but less is known about dynamical correlations. The 
non-local nature of string operator is such that we can not even evaluate simply the lower and upper threshold for 
S and S"' . By applying Sj' to the vacuum, the exponential factor of it will allow creation of any number of 
particle hole excitations 

^+|0>= exp(z.^ctc,)ct|0>= E^^^^%^4lO> 

j<i n 

So it is obvious that potentially S ^{uj,q) could be non-zero for fixed q and any tu. However, in general, we expect 
an infinite number of progressively weaker singularities, extending down to zero energy, similar to the case of the 
longitudinal structure function in a non-zero magnetic fieldii After reviewing the standard bosonization results in 
sub-section A, we discuss effects of irrelevant operators arising from band curvature in sub-section B. 



A. Bosonization approach 

As mentioned before, calculating the dynamical correlations for xxz model is very difficult, and is one of the most 
studied problems of one dimensional spin chains. By using field theory methods known as bosonizationii^i^ we could 
get some information about low-energy effective description of these correlations. Here we review the results of 
bosonization for the transverse spectral function of the xxz Hamiltonian Eq. (jl.ip . We saw that the xxz model 
is equivalent to an interacting fermion model by using the Jordan- Wigner transformation. Now by bosonizing the 
fermion operators we obtain the bosonic representation of spin operators. In the low-energy effective theory, we only 
include excitation around the two Fermi points so fermionic operators can be written in the following form 

c, « ^fl(x)e*'=^^' + iPL{x)e-'''^^ (2.7) 

where ^pn and ipL are slowly varying fields. Now by using Abelian bosonization^! we could write the fields as follow 

4'r,l{^) - -^e-''^-^^.-^^) (2.8) 
V27re 

where e w kp^ is a short-distance cut off, and (^r^l arc the right and left components of bosonic fields (f> and its dual 

field e 

^ (t)L + (t)R 

^ — ;= — 

V2 
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In terms of these bosonic field the effective bosonic description of xxz Hamiltonian is 

1 



n 



LL 



Kid, 



K 



(2.9) 



where v and K are spin velocity and Luttinger parameter, respectively. These parameters are known exactly, for zero 
magnetic field, from Bethe ansatz calculations;^ 



TT \/l — 

v = J-- -, /-iT = [2 - 2arccos(A)/7rl"^ 

2 arccosA ' ^ \ ii \ 



(2.10) 



if = 1 for A = and A' = 1/2 for A 1. For finite magnetic field there is no closed expression for them, but they can 
be evaluated numerically by using Bethe ansatz. The above Hamiltonian is equivalent to the free boson Hamiltonian 
with canonical transformation of bosonic fields as 

(2.11) 



= B = 6/y/K. 

The Hamiltonian then becomes the standard free boson hamiltonian: 



(2.12) 



Now by knowing the Hamiltonian in terms of bosonic field, let's look at the form of spin operators in terms of them. 
This is quite straightforward for 5*1 : 



K . 



S^. « TO+ : -^It/jl : + : f/'Jj'/'i?^ : +(V'iV'i?. + h.c.) = m+ J —dx(t) + C cos \/4t:K(P + (27rm + n)x 



(2.13) 



where : . . . : denotes normal ordering and m =< 5*^ > is the magnetization, which is related to the Fermi wave- vector 
by the exact relation: 



2kp = 27rm + n 

For weak fields, 

m — > Kh/i'Kv). 

is a non- universal constant. To obtain the low energy representation of Sj' from Eq. (|2.ip we use: 

We also approximate the exponential of the Jordan- Wigner string operators using: 

in cjci w ITT / dy{kp + \/ K/irdycj)) = constant -f inkpj + iV TrK(j){j). 



(2.14) 
(2.15) 

(2.16) 
(2.17) 



Following the standard bosonization approach^ we have ignored the oscillating term in CjCj in the exponential of 
the Jordan- Wigner string operator, but we will consider it in the next sub-section. Note that exp[i ^^-^^ cjcj] is 
Hermitian, taking eigenvalues ±1. On the other hand, the exponential of the continuum limit operator in Eq. (|2.17p 
is not Hermitian. To deal with this problem the standard approach^ is to instead take the continuum limit of 
cos[7rX;i<j4c/]: 



cos[7r^c|Q] oc e^^kirj^iV^^U) ^ ^^^^ 

l<3 

Substituting these low energy limit formulas into the second of Eq. (|2.ip gives 

C{-iy + C- cos(27rmx + ^/A-KK(j)) 



(2.18) 



SJ oc (-l)Je-V-/^« 



(2.19) 
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(a) S+- 



(b) s- 



FIG. 1: Singular points of zero temperature transverse spin spectral function of xxz model, predicted by 
bosonization. The solid lines indicate diverging singularities and dashed line the vanishing singularities. 



where C and C_ are two other non-universal constants. It can be seen from Eqs. (|2.13p and (|2.19p that the applied 
magnetic field induces a shift in momenta of both transverse and longitudinal spin operators, but in different ways. 
For the longitudinal operator, it shifts only the staggered part but for transverse spin it shifts only the uniform part. 
This shift is small for weak fields where the approximation of ignoring band curvature is valid. 

Let us focus on transverse structure function near g w for weak fields. At zero temperature we have 



^—iHx 



{vt-x- «e)2+2';(wt + x- ie)2'; 

where x — j, H and r/ are given by 

H = 27rm 



(2.20) 



ri^^ '- 2.21 

' 8K ^ ' 

obeying 77 < 1/8 for |A| < 1, and e is a positive quantity of order the lattice spacing. By taking the Fourier transform, 
the spectral function is 

b {q.Lo) ca 0{uj — v\q + Hy 



{u-v\q-H\) 



{io - viq + H)y-'^v 

{LO-v{q-H)Y+^'^ 
{u} + v{q - H)y-^-n ' 



Note that the first term has a diverging threshold at a; = v{q + H) ior q + H > and a vanishing threshold at 
w = —v{q + H) for q + H < 0. The second term is the parity transform {q — >■ —q) of the first. Diverging and vanishing 
thresholds are indicated by solid and dotted lines in Fig. [1^. Note that, two diverging thresholds occur in 5*^ (q, uj) 
for l^l < H which cross each other at g = but that there is only one diverging thresholds for |g| > H. S ^ is 
obtained from by the transformation H — > —H. (Recall that we are assuming H > 0.) 



S (q, 0.) ce{u-v\q-H\) ^^_ + ~ v\q + H\) ^ ^ (2.22) 

Its diverging and vanishing thresholds are shown in Fig. |lb| Note that, for S ^, no diverging thresholds occur for 
\q\ < H and a single diverging threshold occurs for \q\ > H . 
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FIG. 2: Zero temperature transverse spectral function {uj, q) predicted by bosonization for fixed q. a) shows 
transverse spectral function for zero magnetic field and b) is for non-zero magnetic field with q < H . 



For ESR applications we will be especially interested in the case A slightly less than 1 and small H corresponding 
to K slightly greater than 1/2 and thus 77 ^ 1. Then it is important to note that the 77 dependence of the constant 
factor in and S is^^ oc sin'^(27r77)r(— 1 — 277)r(l — 277) which vanishes linearly with 77. Here F is Euler's Gamma 
function. To study the line shape at _ff = and 77 0, we take into account that this expression for S{q, w) is only 
valid for a finite range of w, 7;|<7| < w < A for an upper cut off A, of order J or less. We then use the fact that 

Ih, /\'^^^^-/^) ^hmA^"^!. (2.23) 

Therefore we have 

2r]0{u-vq) 
um ,.. „ = d(uj — vq) 

and thus the term with a diverging threshold approaches 

S{q,Uj)(Xv\q\5{L0~vq). {rj^Q). (2.24) 

For a fixed momentum q and small 7/, S{q, u) as function of lo is depicted in Fig. [J. Here we show only one term 
with a diverging threshold. It is zero for frequencies such that ui < vq, and it has a local minima at point w* given by 



'-^vq^^ (2.25) 

So for 77 « 0.1 we get w* « 3vq. We should also be careful about the cases of very small anisotropy; from Eq (|2.25|) 
we see that for fixed momentum as anisotropy gets smaller and smaller oj* becomes larger and larger, so it seems 
that we are getting out of the region where bosonization works. The results of bosonization are reliable below some 
cutoff A; then the consistency relation uj* < A will gives us a restriction on momentum q such that we must have 
vq < y^A. 

In fact, as shown above 5"^ and S ^ are a sum of two terms each with a separate threshold for all q 0. 
Depending on which spectral function we look at and the value of q, these thresholds can be both diverging, both 
vanishing or one diverging, one vanishing. The various shapes of S{q,uj) are sketched in Fig. [3]. In the special case 
q = 0, there is a single term of diverging threshold type. 

For the transverse Green's function at q sa tt we have 

< SH^,t)S-iO) >oc< S-i.,t)SHO,0) >oc ,)v^l^,„,,)r/.. (2-26) 



CO, CO, 



CO, CO, 



(a) \q\ < H for 5+" 



(b) H < \q\ for 5+" 





(c) |g| < H for 5" + 



(d) H < \q\ for 5" 



FIG. 3: Various shapes of singularities of 5"+ and S ^ for different range of momentum q, with \q\, _ff <C 1, 

predicted by bosonization. 



where x ^ j. By taking the Fourier transform we obtain the transverse spectral function near g = tt: 

9{u! — v\q — tt\) 



S {q,(^) ~ S '^{q,uj) (X 



- 7r^2^l-l/4^!: 



(2.27) 



So we see that the singularity exponent for staggered part is 1 — l/4:K, which is different than the exponent for 
uniform part. Now there is a single diverging threshold at co ~ w|g — 7r| for either sign of g — tt for both S'^ and S ^ . 
It seems natural to assume that the diverging thresholds of S'^ and S starting at q — ±H can be interpolated to 
the diverging thresholds terminating at g = tt. Such an assumption goes beyond the standard bosonization approach 
which is restricted to q near and tt (and to small H) but we shall see in the next section, using X-ray edge methods, 
that this interpolation is correct. However we will also find that the exponents of the diverging thresholds are predicted 
by standard bosonization; 1 — 2rj near \q\ = H and 1 — 1/{4:K) near q = n are both incorrect, as are the exponents 
of the vanishing thresholds. 



B. Effect of irrelevant band curvature operators 



In previous sub-section we reviewed the prediction of bosonization for the singularity exponent of transverse spectral 
functions. In this section by including the effect of band curvature operators we show that the predictions of naive 
bosonization for singularity exponents are not reliable. 
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In bosonization approach we treat the interactions exactly but we linearize the dispersion around Fermi points and 
neglect the effects of higher order band curvature terms. By power counting these terms are irrelevant in low energies 
and renormalize to zero, but as discussed in the effect of these operators is important near the singular thresholds 
for longitudinal spectral function. We will show that the same argument works for transverse spectral function. In 
this section we look at the effect of these terms to lowest order for h ^ 0. 

As shown in [13 and by including the effect of band curvature corrections the Hamiltonian becomes 

H = Hll+SH (2.28) 
where Ull is Luttinger Liquid Hamiltonian and 61-1 is given by 



+ T]+\{d^.^Lfd.^R - {d.^Rfd^^LW (2.29) 



Where to first order in A we have 



^ ( 2A . , 

?7_ « — 1 H sm Kp 

ra \ TT 

3^ • I. 

77+ w sm kp 

■Km 

Where m = (J cos kp)~^ is the effective mass of Fermi excitations. For weak interaction we can neglect r/^ and only 
include the effect of rj- term in Eq. p.29p . 

Now let us evaluate the transverse spectral function using perturbation theory in SH. By ignoring terms proportional 
to 77_|_ which mixes right and left operators, in general we are looking for the following kind of imaginary time correlation 
functions 

G<.i/(a;, r) = G'_R(a;, r)G'L(x, r) 

Gr{x,t) = < e*N/2¥i:vR(a;,r)g-»\/2¥i;<PH(0,0) ^ 

Gl{x, t) = < e»^/2¥Pv>I,(a^r)g-^^/2¥PvI,(0,0) > (2.30) 

Where h' and v could be written explicitly in terms of Luttinger parameter, K, but for the following discussion we do 
not need their explicit form. The first order correction from perturbation Eq. (|2.29p modifies the correlation function 

to 



(2.31) 



To evaluate the correlation function Eq. (|2.3ip . we focus on corrections to Gji{x,T); calculations for Gl{x,t) are 
exactly the same. 
So we have 

GR{x,T)=Gf{x,T) 



'Jh:pl j dzdr' < e'v^^«(-.-)(a,(p^(z,r'))=^e-^^^«(°'") >o (x,t) (2.32) 



In diagrammatic way the non-zero contributions are depicted in Fig. [3] 
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(0,0) 




(0,0) 



FIG. 4: The relevant Feynman graphs which contribute to correlation in first order of perturbation. 



So the non-zero correction is given by 

6Gr (X J dzdr' < e*^^'^^(^^^)(a,(pflX^,^'))^e-*^^'^«(°^°) > 

I ^ — ^ r?, ' rr?, ' 



/ dzdr' > 3n m m - 1)— ^ ^ 

7 ^ ' nlnil 



<^k(x,t)"-V«(0,0)™-2> 



X < ipR{x,T)d:,ipR{z,T') >< d^(fR{z,T')ipR{0,0) 



(2.33) 



In going from second line to third line, we have used Wick's theorem, and the factor in m{m — 1) comes from the all 
possible number of contraction of fields. Thus we have 

6Gr cx -3^V(27^^.)3Gg^(x,r) / dzdr' < (^^(.t, t)9,(^^(z, r') >< 9,<^fl(z, r')¥>ii(0, 0) >^ (2.34) 



Now by using the fact that 



< dx(pR,Lix,T)ipRx{0,0) >= 



1 1 

27r VT =p ix 



(2.35) 



we can write Eq. p.34p as 



6Gr oc 3iJ(2ni'fG^°\x,T) [ dzdr'— ^ 

I 111 n-' — -T- 1 — 



v(t' — r) — i{z — x) [vt' — izY 

= -i2.zy(2^--L_ / - - 

(uT — ixY 
= -127riv/(27rt^ ^ 



[vT — ix)'^ 



By replacing vt {{vt — ix) + {vt + ix))/2 we can write SGr as 



(2.36) 



Xri ^ S'^V (27rt/)3 f 1 . vT+ix \ 

OUR (X - l^(„r-ix)i + - + J 

With exactly the same calculations we will get results for left moving fields as 

67r-i-y/(27rP)3 f ^ , vr-ix \ 

"^LO^ 1, y{vT+ix)^+^ {vT+ixy^+" J 
Now by plugging all these results into Eq. (|2.3ip . the final form of the correlation function to first order is given by 



(2.37) 



(2.38) 



G^i^{x, t) 



-{1- 



^ VT + ix 

{vt ~ ixY {vT + ixY V \VT ~ ix {vt ^ ixY 

47r^77_-\/F^ / 1 VT — ix 



VT + ix {vt + ix)'' 



} 



(2.39) 
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Now by taking the Fourier transform of Eq. (|2.39p and continuing to real frequencies, we have 

V \ LJ + vq J 

+ {LO + vq)[l + ^ } 2.40 

V \ uj — vq J 

Where ry- w 1/m. It is easily seen from the last term of Eq. ()2.40p that, as uj approaches vq, the pcrturbative 
corrections blow up like q^ /m{uj — vq). This is exactly the reason that bosonization fails near the threshold. We 
also see that for \uj — vq\ ^ q^ /2m, perturbative corrections become small and we get the naive bosonization results, 
so we showed that irrelevant operators potentially will change the singularity exponent of correlation functions near 
the threshold but their effect is negligible away from the threshold so they would not change the qualitative shape 
of correlation function found by bosonization. In the following section by using the X-Ray edge method we find the 
singularity exponents of transverse spectral functions. 



III. X-RAY EDGE METHOD 



In sub-section III Al by use of standard bosonization we found that the transverse spectral functions 5*^ (w, q), for q 
near has a diverging singularity with exponent 1 — 277. In sub-section III Bl we argued that band curvature operators 
would change the result of bosonization for singularity exponents. In this section we will explore this question by 
extending bosonization using X-ray edge methods introduced in d-d]. We find that by use of these methods, transverse 
spectral functions have different critical exponents at singular energies than predicted by standard bosonization. We 
also point out the existence of large numbers of sub-leading singularities with vanishing intensities, similar to the 
ones found by standard bosonization. We use the notation and results of and we skip the details of derivations; 
interested readers should see and 01 for most detailed calculations. In this approach we try to evaluate transverse 
spectral functions for fixed momentum q by including the effect of a single high-energy particle or hole excitation. To 
do this we need to find relevant momenta of this excitation which contributes to the spectral functions at momentum 
q. In this section we find the effective Hamiltonian for these excitations, and also those relevant momenta. Suppose 
that the momentum of this particle or hole excitation is k. (We will eventually use the notation fcp for a particle and 
kh for a hole.) In X-Ray edge method we are interested in the high energy excitations near this momentum k and 
also low-energy excitations around the fermi points. Thus wc can write the fermion operator in the following form 

Cj- « Vi?e''''^^' + i^Le-"-^"' + d e'^^ . (3.1) 

Where ipRjipL and d vary slowly on lattice scale. Then by linearizing the dispersion around the fermi points and 
bosonizing low energy fermions, and also linearizing the dispersion around high-energy particle or hole excitation, we 
have 

v 1 

U = d^ [e- iudx) d + - [(S^j^^l)^ -t- {OxIPr)^] + —7== {^Ldx^L - undxipR,) d'^d 
^ V 27riv 

Here the chiral fields Lp^ and lpr are the transformed ones defined by: 

Q ^ + ^R 

This Hamiltonian is described in [ [26l - l28l | for Luttinger liquid coupled to an impurity. The parameters of above 
hamiltonian are as follows; e is the energy of the high energy particle or hole and, for A = 0, it is given by e = 
— 2 J(cos k — cos kp). u is the velocity of the heavy particle or hole and equals J sin fc at A = 0. v and K arc the boson 
velocity and Luttinger parameter, respectively, v may be regarded as the Fermi velocity of the interacting fermion 
model. It is the only velocity appearing in the standard bosonization approach and plays the role of the "velocity of 
light" in the effective Lorentz invariant field theory. The velocity parameter u, describing the high energy particle 
or hole is an important new parameter in the X-ray edge approach. Finally kr^l are the couplings between high 
energy fermion and bosonic fields and to first order in A are given by kr^l — 2A[1 — cos{kp =F k)]. These coupling 
could be evaluated by Bethe ansatz calculations^ for any A, H and k. They are important for finding the singularity 
exponents. Note that, in general, u, e and kr^l all depend on k as well as A and h. 
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The Hamiltonian p. 21) looks complicated as it contains interactions between fermions and bosons. We can eliminate 
the interacting part of p.2p by doing a unitary transformation given by 



U = exp 



dx 



V2^ 

Effect of unitary transformation on bosonic and fermionic fields is as follows 



(3.2) 



d = d exp 

VR,L = 'pR,L ± 



1R,L 



=d\x)d{x) 



V2t:K 

1R,L 



N{x) 



where N{x) is defined by: 



N{x) 



4V2^ 



sgn(a; - y)d'^ {y)d{y)dy. 



It is easy to see that unitary transformation leaves d'^{x)d{x) invariant. With 7ij.L given by 

I^RX 



1R,L 



V ^ u 



(3.3) 



(3.4) 



wc can decouple fermionic fields from the bosonic ones. Having done the unitary transformation the Hamiltonian will 
look like 



H = - [{d^^Lf + [d.^Rf] + (e - ma,) J - 



(3.5) 



Where • • • means higher dimension irrelevant interactions that will be produced by doing unitary transformation and 
which we ignore. 

We now consider the transverse Green's function: 

s~+ =< S-{t)S+iO) > . 

By doing Jordan- Wigner transformation, we have 

S~+{t) = e'""^ < Cj{t) cos[t:N J (t)] cos[7riVo(0)]cS(0) > 

where Cj(t) is approximated as in Eq. p.l|) and 

To obtain the transverse spectral function, 5* ^{q,u;) at a wave- vector q far from the low energy regions, ±iJ, tt, the 
term that we are interested in is 



S~+it) = e"^'e^'=''J' < d(j,t)cos[7r7Vj(i)]cos[7riVo(0)]dt(0,0) > . 



(3.6) 



We see that d must be chosen to be a particle operator and we have consequently labelled its momentum kp. Note 
that we have written the Jordan- Wigner string operator in manifestly Hermitian cos form, as in Sec. IIIA. Now we 
decompose Nj{t), into c- number, non-oscillatory and oscillatory parts 

Njit) = —j + n{j,t)+m{j,t) 

TT 



h{x, t) 
m{x, t) 



dy : ijl{y,t)i,'R{y,t) : + : tPl{y,t)ijL{y,t) : +d^y,t)d{y,t) 
^'^^{y, t)d{y, t)e'('^^-^-)y + ^{{y, t)d{y, t)e'('^-+'^-)y + ^(y, t)i^R{y, t)e''^-y + h. 



i{kp — kp) 



^ {x\t)d{x',t)e^^''^-'^-^^ + ^^^^[(x',t)d(x',t)e^(^--+'^>)- + -^i;l{x' ,t)ijR{x' ,t)e'^'^-^' + 

i[Kp + kp) Zikp 
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where e — O"*" , x' = x — e and in the third line we used the fact that both ij) and d are slowly varying fields and most 
of the contribution of the integral conies from limiting point x — e. At this point we will set the rapidly oscillating 
term, TO(a;, t), to zero. This will give the dominant divergent singularity in the transverse spectral function. By Taylor 
expanding in powers of m(x,t) we obtain various vanishing singularities as well as unimportant renormalizations of 
the amplitude of the divergent singularity, as we discuss in Sec IV. We may then decompose h(x, t) into its commuting 
high energy and Fermi surface part. 

h{x, t) = n{x, t) + nci{x, t) 

n{x,t)= / : Tpl^{y,t)TjjR{y,t) : + : ipl{y,t)ijLiy,t) : 

J OO 

nx~~e 

ndix,t)= / d^y,t)d{y,t)dy. 

J oo 

Because all the d operators in h{x, t) are at points y < x, we have 

[n{x,t),d{x,t)] = 

and thus we may drop the Ud terms leaving: 

S-+{t) (X e^f'^'f+'^^J' < coslkpj + 7Tn{j, t)]d{x, t)d^(0, 0) cos[7rn(0, 0)] > 

Following Eq. (|2.17p this becomes: 

S'j^+(t) cx e't'^p+'^^J' < cos[fcFj + V^(t>ij,t)]d{j,t)d'' (0,0) cos[\/^(l){0,0)] > . 

Note that we have treated the Jordan- Wigner string operator in precisely the same approximation as in the standard 
bosonization approach. We now make the unitary transformation of Eq. (j3.3p so that the fermion and bosons are 
decoupled. Noting that N{x,t) annihilates the vacuum this leaves: 

S-+ oc < d>,i)cos[fcFa; + V^^(x,i)]e^f^"^"^"'*^+^"^""^"'*" 

ev*s[T«'^«("^°)+^-'^"-(°'")lcos[V^0(O,O)]Jt(O,O) > e*(^>+-)- 

Separating the fermionic and bosonic factors, this becomes: 

Sr+{t) = 5W-+(a;, t) < d{x, t)d{Q)^ > (3.7) 



where x = j, 



and 



= e'''''+''''-''^^Tix,t)+e'^''+''''+''^^''I+ix,t) (3.. 



/+(a; t) = ^g-i^/'^'^'^RVR(x)+i^y2^Tly+tpL{x)^^^/2^^ly+ipR{0)-^^y2^^ly+tpL{0) ^ 



^ £ + ivt ~ ix J \ e + ™t + ix , 

1~ [x t) = V 27ri/^i^R(a;)-i-y/ ix) g-^V 27ri/^<PH(0)+i-^ 2171^^^ ipi,(0) ^ 

e 



e + ivt ^ ix J \ e + ivt + ix 
Where e is of order of the lattice spacing and ^ are defined as follow 



(3.9) 
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At zero magnetic field we have ^r^l/t^ = 1 ^ independent of momentum, so tlie results will simplify to 



4 = —{l-K±Kf 

4^^{l-KTKf. (3.11) 



Having done the unitary transformation, the d fields act as free particle so we have 

< J(x,i)Jt(0) >« e-'^^"^* I ^e''^^^-"*) « e-'^^'-^'Six^ut). (3.12) 

We can now turn to the question of how the momentum of the high energy particle, fcp, should be chosen to study 
S ^(g, Lu) near the threshold for arbitrary q. From Eq. p.Sp we sec that a high energy particle of momentum kp gives 
terms in the transverse Green's function oscillating at wave- vectors t: + kp — kp and tt + kp + kp. Thus we see that 
there may actually be two choices for kp which will give a contribution to the transverse Green's function oscillating 
at a specified wave-vector q: 

k^ = q + TrTkF- (3.13) 

In general, both must be considered in calculating the singularity behaviour of G ^{q,Lu). However, the /cp's are 
restricted by the requirement that they are allowed particle momenta, 

kp < kp < 2tt - kp, (mod 2tt). (3.14) 

(Recall that kp — tt/2 + ttto > tt/2.) Thus we see that the high energy particle of momentum fc+ contributes to 
S ^{q,Lu) for q in the range [2kp — TT, n] and the particle of momentum contributes for q in the range [— tt, TT — 2kp]. 
For a given q there is at most one possible high energy particle momentum contributing to 5* ^{q, co). Having identified 
the appropriate high energy particle momentum we may now complete the calculation by Fourier transforming Eq. 
([XT)) using Eqs. (jX^ - ipi^ . The resuh is: 

p f.il^-<K')^*S(x - ut) 

5'"+((j,Cj)(x dxdt — i {2kp - TT < q < tt) 

J [vt — X — leyR {vt + X — ieY^ 

(X dxdt ^ {-Ti<q<ii -2kp). (3.15) 



[vt ~ X — ie^R {vt + X — ieY^ 

The x- integrals may be done trivially using the (5-functions. 

r gi["-«(fcj)]t 

S ^((7,cj)cx / dt T ^, (2kp — IT < q < -k) 

J [{v ~ u)t - ieYn[{v + u)t - leyi 

g''["-<!('=p )]* 



dt {-t: < q <TT ~2kp). (3.16) 

J [{v - u)t - ie^R [{v + u)t - iefL 

The ^-integrals can now be done by contour methods. Note that, if t; > u, they are only non-zero for ui > i{kp) 
and uj > e{k~) respectively, corresponding to a lower threshold. Using the free particle cosine dispersion relation 

V > u is always satisfied for particle excitations. While this dispersion relation is known to be exact, apart from an 
overall factor, including interactions for ft. = it is in general modified. We might expect that v > u remains true for 
particles, at least for small enough h. However, see below. Assuming this, we obtain: 

S ^{q,oj)(x- (3.17) 

[lj - u}L(q)r^''' 

where the singular energies are given by: 

iOLiq) = eikp), (2kp - TT < q < it) 

= e(fc~), {-TT<q<TT-2kp). (3.18) 

The critical exponents are given by 

ti{q) = 1 - v^{k+) - vl{k+), {2kp -TT <q<n) 

^l~v^{k-)-uJ^{k-), {-n < q < TT - 2kp) (3.19) 
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with ]^{k) given by Eqs. p.lOp and (|3.11|) . Note that the same phase shift parameters Jr/l determine both 
longitudinal and transverse spectral functions; however we need to know them at the momentum of the high energy 
particle or hole which is not the same for longitudinal and transverse spectral functions, for given q. In general, 
the phase shift parameters depend on momentum as well as field, becoming momentum independent at /i = 0. It 
follows from parity that ^^(fc) = Hft{—k). Since, for g > 0, kp{—q) = —k^{q), (mod 27r) g) = ^^{<l) and 

^Ri-(l) = ^til) and hence ^i{-q) = ^i(q). 

It is interesting to compare both the singular energies and exponents to those predicted by standard bosonization as 
q approaches the zero energy points q w ±H = ±{2kF — tt) and q ^ n. Near these zero energy points we may linearize 
the e(fc) giving: u>Li<l) ~ ''^Iq-F H\ and v\q — 7r|, precisely the singular energies predicted by standard bosonization. 
This X-ray edge calculation also confirms the conjecture made in Sec. II that the diverging singular energies at the 
different low energy momenta are smoothly connected. On the other hand, the exponents appear to disagree with 
the standard bosonization results for all A and h, a similar observation to the one in Ref. This can be seen, for 
example, by considering the limit A — 0. In this case 7^ ^ — J' 0, K ^ 1 so j/^^ — > 1/4 and — >■ 1/2, for all g and h. 
On the other hand the standard bosonization result from Eq. (|2.2ip and (|2.22p is 77 1/8 and hence /i ^ 3/4. We 
can also compare the zero field predictions for general A. From Eq. (|3.1ip 

lj = 2-l/{2K)~K, {h = Oyq). (3.20) 

On the other hand, standard bosonization predicts yu = l — 277 = 2— 1/{4K) — K for g « and = 1 — 1/{4:K) for 
g sa TT. We expect that standard bosonization fails to predict critical exponents correctly for the transverse spectral 
function, as discussed in sub-section IIB. Given this situation, it is useful to check the SU(2) symmetric case, h ~ 0, 
A = 1. In this case X-ray edge methods predict, from Eq. (|3.20p = = 0, = vj^ = 1/2, /i = 1/2 for the 
transverse spectral function independent of g. The same exponents were found earlier— for the longitudinal spectral 
functioni^ (In this case, they agree with standard bosonization near g = tt but not near g = 0.) 

By doing similar calculations we could find S'"* , which is different than S ^ for h ^ Q. Following the same 
procedure, Eq. (|3.7p is replaced by: 

S+- (X < J^(:E,i)cos[/fcfX + V^0(x,<)]e^'^'''^'*^"'*^+^^''"^*"'*^' 

^-^[iHVH{m+iLVL(m] cosV^0(o,o)]J(o,o) > e*(-'^-'>+-)- 

where d now annihilates a particle in a filled state below the Fermi energy with momentum fc/i, i.e. creates a hole. 
This can again be factorized as: 

S+-{t) = 5(°)+-(a;,t) < d\x,t)d{Q) > (3.21) 

where 

S^°^+-{x,t) ^ e*(^-'=''+'=^)^/-(2;,t) 4-e'('^-''"''-'=^)^/+(2;,t) (3.22) 

and I^{x,t) are the same functions defined in Eq. p.9p . except that kh must lie in a different range, \kh\ < kp. Thus 
defining: 

kl = t: - q^kp (3.23) 

we see that the first term in Eq. p.22p is non-zero for g in the range [tt— 2fc^, tt] while the second is non-zero for q in the 
range [— tt, 2kF — tt]. These are wider ranges than occur in S ^ . (Recall that we assume H > and hence kp > tt/2.) 
In particular, both terms can contribute for |g| < H. Again, as g approaches the zero energy points, ±H and tt the 
singular energies approach those predicted by standard bosonization. Again as anticipated in sub-section IIA, the 
singular energies at these zero energy points can be smoothly connected. Another interesting feature is the shape of 
the singularity. For S ^ the singularity was one-sided, vanishing for lu < u^. This was a consequence of the fact that 
the velocity of the high energy particle always obeys u < v assuming this feature of the non-interacting dispersion 
relation is unchanged by interactions. On the other hand for holes, again using the non-interacting dispersion relation, 
u < V is only obeyed if |fe/i| < tt — kp; there is a range of hole momentum near kp where the high energy hole has a 
higher velocity than the Fermi velocity. \k^\ < n — kp corresponding to < g < 27r — 2kp ^ n — H. Thus, in this 
region the singularity is one-sided, oc 0(w — col)- On the other hand for —H < g < and tt — H < q < tt, where 
u > V the integral in Eq. p.l6p is also non-zero and gives the same critical exponent with a different amplitude for 
Lo < ujl. In this case we see that uj^^q) is not a lower threshold. There is also spectral weight below this frequency. 
The qualitative shape of 5'"' (g, w) for < g < is sketched in Fig. [5] . 



16 



We emphasize that a singularity is one-sided for u < v and two-sided for u > v where v is the Fermi vefocity and u 
is the velocity of the high energy particle or hole. In general, these velocities depend on A and also /i, being strongly 
renormalized by interactions. In [IJ], Fig. [15], it was illustrated that for A = 1 and small non-zero field one of the 
cubic term in the bosonized Hamiltonian, due to band curvature effects, has a coupling constant 77_ < 0. This may 
indicate a reversal of the sign of the effective mass due to (strong) interaction effects, implying a reversal of the sign of 
u — V as the energy of the high energy particle or hole approaches the Fermi energy, corresponding to \q\ — > \2kp — ■k\ 
(or g — >■ tt). Thus in this parameter range, the one-sided and two-sided nature of the 3 singularities in and S ^, 
discussed above, would be reversed. 



FIG. 5: The behaviour of 5*"*" (a;,*;) correlation for |g| < H. There are two different hole excitations which 
contribute to the spectral function; the lower energy hole produces a cusp-like singularity. 

Given that the singular energy, wl('?) is sitting inside a region of non-zero spectral weight, for certain ranges 
of we might ask whether it is reasonable to expect singular behaviour at this energy or whether the infinite 
peak might be broadened and made finite due to some sort of decay process for this high energy excitation. This 
important question also arises for the longitudinal structure function and for the fermion spectral function. It has 
been suggestediii^ that integrability might prevent this excitation from decaying, for some range of momentum, 
even though it is kinematically allowed, leaving the singularity intact. This is true because three-body scattering 
processes are required for it to decay and these are expected not to occur in this integrable model. This seemed to be 
consistent with Density Matrix Renormalization Group results for the fermionic spectral functionii In fact, we should 
also consider processes in which the heavy particle or hole decays by producing 3 other high energy quasi-particles, a 
2-body process which is expected to be present even in this integrable model. Using the — cos k dispersion relation, 
this is kinematically allowed for high energy holea^ in the region u < v, corresponding toO< \q\ < n — H (the upper 
branch for < |g| < H) but not allowed for high energy particles for any g.— Putting these observations together, 
we expect the sharp one-sided singularity of S ^(q,uj) to be present for H < \q\ < t: and the 2-sidcd singularity of 
{q,uj) for < \q\ < H and tt — H < \q\ < n to he present. The one-sided singularity of for if < \q\ < tt — H 
should be broadened by higher order interaction effects^ not taken into account in this treatment. See Fig. [B]. Again 
we emphasize that the precise region of q over which singularities are broadened depends on the dispersion relation, 
which is modified by interactions; here we have just stated it using the — cosfc dispersion relation, valid at small A. 

So far, we have set the rapidly oscillating operator m{x,t), defined in Eq. (|3.7p to zero. The effects of including it 
are discussed in the following section. It basically leads to additional terms in the transverse spectral function which 
have singularities at different energies, including lower ones. However, these singularities are of vanishing type, with 
exponent < 0, dashed lines of Fig. [S]. The relatively simple approach we have taken here is just sufficient to give 



S i 
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(a) S+- (b) S-+ 

FIG. 6: Singular points of zero temperature transverse spectral functions of the xxz model, predicted by X-ray edge 
method. The solid lines indicate diverging singularities and dashed line the vanishing singularities. The grey thick 
line indicates the range of momentum that interactions should broaden the diverging singularity. In Fig. [S^ and 
|6b) the lower dashed line is given by particle excitations and upper one is given by hole excitation. 



the diverging singular terms. 

The situation is considerably simpler at zero field, /i = 0. In this case, the free dispersion relation is known to be 
exact, apart from an overall change of amplitude, 2t — > v. Thus the condition u < v is always satisfied so S ^ = 
has only one single-sided singularity at ujl = v sin q with v given in Eq. (j2.10p and critical exponent given by Eq. 
(|3.20[) . In this case, no decay processes are kinematically allowed and no additional singularities occur, since the single 
hole or particle has the lowest possible energy for given wave- vector. 



IV. SUB-DOMINANT SINGULARITIES 



In Sec. II and III we ignored rapidly oscillating terms, m{x) of Eq. (j3.6[) in the Jordan- Wigner string operator 
in calculating the transverse structure function and also the effects of Umkalpp term which oscillates as e''*'^^^. We 
consider the effect of including these terms here. Let us begin with the term: 



Zikp 



(4.1) 



Actually, this term represents a correction to standard bosonization, even without using X-ray edge methods, so we 
consider its effects there. To make things as simple as possible we also consider zero field, kp = 7r/2. Then, after 
bosonizing m, the standard bosonized expression for in Eq. (j2.19p is modified to: 



57 oc e-'V^e(-) 



C{-iy + cos(V47rAX.T)) exp[(-l)J'iC" sitl{^J ^-n K (l){x))] 



(4.2) 



for a non- universal constant C . We now Taylor expand the exponential and use double angle formulas. Wc see that 
the staggered and uniform parts of SJ have a series in increasingly irrelevant operators: 



o„ — e V / ' ' a2„+i( 



i2n\/iTTK4> 



,i{2n+l)V4:TTKct> 



(4.3) 



The effect of including the e*™ factor is simply to renormalize the coefficient of the leading operator in and Sg 
together with producing the irrelevant corrections. Now consider non-zero field. Eq. (|4.2p gets replaced by: 



Sy cx e-'V^»(^) \c{-iy - C- cos{{2kp - 7r)j + VM(/)(x))j exp[^C"e2*'=^^'+*^'^'^(^) + h.c.]. (4.4) 
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We again get a series of irrelevant operators but now all at different wave-vectors: 

C- _ -iJ IT I KB(x) i(2nfcF+7r)j+in\/4irK0 
— ^ ^ ^ tint- 



(4.5) 



It is interesting to note that this expansion contains precisely the same terms as the one derived by Haldane^^ for a 
boson annihilation operator in a Luttinger liquid. We also see why the replacement of the exponential Jordan- Wigner 
string operator by a cosine form, its Hermitian part, is not really necessary. Keeping the complete expansion in Eq. 
(14. 5p . we get the same set of operators either way. 

Now consider the effect of the term in Eq. (|4.1[) in the X-ray edge approach. After the unitary transformation of 
Eq. p.3p the term in S~ linear in the d operator, with momentum fc, is: 



S, (X de ' 



,i[{2n+l)kF+T+k]j+i(2n+l)\/iTK4> 



(4.6) 



This expansion is similar to the one derived by Haldano^ for a fcrmion annihilation operator. The momentum q at 
which the n*'* term contributes to the spectral function is: 



{2n+l)kF + TT + k 



(4.7) 



where the momentum k must correspond to that of a high energy hole, < kp in calculating or to that of a 
high energy particle, kp < \k\ < tt, in calculating 5* '". Note that by ignoring m{x) in Sec. JIB we only considered 
the n = and n — —1 terms in the sum of Eq. (|4.6|) . The n*'' term in the expansion of SJ in Eq. (|4.6|) leads to a 
singular term in the transverse spectral function. 



S{q,uj) cx 



1 



\LJ-Uj{q)\t^ 

at the energy uj{q) given by the energy of the corresponding particle or hole: 

a;„(g) = ±e[q - {2n + l)kp - tt]. 



(4.8) 



(4.9) 



We see that, for general kp, these thresholds can occur at arbitrarily low energy. (For rational kp there is a finite 
number of them and there a lowest one at a non-zero energy.) The corresponding critical exponent is given by: 



,(") 



(4.10) 



where z/j^"]? are the left and right scaling dimension of the n*'* operator in Eq. (j4.6|) . These obey: 



in) 



in) 



IB. 



(2n+ 1)VX + 



{2n + l)y/K 



IR - IL 



2 



= - {2n + 1)VA' 



2'K^JK 



(4.11) 



Thus the two largest values of n occur for n = and —1 for [7^^ — ^]^\/ {2ttK) < 1. At zero field, 7^ = 7^^ so this 
condition is satisfied. Also at small A, ''yL,R are 0(A) so the condition is again satisfied. It should remain satisfied 
for a large range of field and A quite possibly including the entire Luttinger liquid regime, but without determining 
the 7L,fl explicitly we can't determine this range. It certainly includes weak fields h <^ J relevant to most ESR 
experiments. For zero field we have explicitly: 



,(") 



K , 



2K 



(4.12) 



so we see that the exponents for sub-dominant singularities (n ^ 0, —1) obey /i < —3/2 for all K in the Luttinger 
liquid regime ii:>l/2. AtA = 0, 



(ti) (n) 



-(2n+l)' 



(4.13) 
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and fi < —7/2 for all sub-dominant singularities. We expect that ^ < for all sub-dominant singularities a wide 
range of field and A including the weak field regime. So the approximation of dropping m{x) made in sub-section IIB 
appears quite generally valid. 

However, we must also consider the other terms in to cx i^jid, i^\d. These give contributions to the transverse 
spectral functions proportional to matrix elements in the d space containing more operators. Consider, for example, 
the case of a high energy particle. Then, to first order in these operators we either obtain d(— e)(i^(0)|0 >= or else 
a matrix element: 

< 0\d{x, t)d{x - e, t)d) {-f)d\{))\Q > . (4.14) 

In fact, this is also zero as e — > as follows from Wick's theorem and translation invariance: 

< 0|J(a;, t)d{x - e, i)d^(-e)d^(0)|0 >=< 0\d{x, t)d'' {0, 0)|0 >^ - < Q\d{x - e, i)J^(0, 0)|0 >'^^ 0. (4.15) 

Now consider the higher order expansion in the terms in to proportional to d and d^ . Even orders in the expansion 
give 

[J(-e)Jt(-e)]»Jt(o)|o >oc d^0)\0 > . (4.16) 

The factor in the Green's function involving Fermi surface excitations has S~ dressed by n ipL,R operators and n 

operators, all at the same point. The {tp^ipR)^ and (V'flV'L)" terms just give contributions the same as Eq. (|4.6p . 
modifying the coefficients. Other products give higher dimension operators using: 

oc A + Ba,0fl-f Ca,(/)L (4.17) 

et cetera. Odd terms in the expansion in the terms in m, which are proportional to J and d\ give zero as e — >■ 0, since 
they are proportional to the same matrix element, Eq. (|4.14[) . 

Till now we have considered the effect of particle excitation on S '". In general high-energy hole excitations, also 
could contribute to S as well as particle excitations to S'^ . Now we will find the effect of hole excitations on 
S ^ . Hole excitations only give us higher order corrections to S ^ , which have sub-dominant singularities, but to 
complete our discussion wc find the vanishing singularity exponent of hole excitation to S similar argument holds 
for particle excitation effects on . In addition to the Eq. p.6p . there is another term which contributes to S 
of the following form 

Sj+{t) = ^i^n < ^^(j- t)cos[^iV,(t)] cos[7riVo(0)]Vlj(0,0) > . (4.18) 

by similar argument as discussed in section IIIII wc could do the canonical and unitary transformations to decouple 
bosonic fields from high-energy excitations. Now the zero order term gives us the naive bosonization results; the 
interesting contribution comes from the first order expansion of to oc tp]^dh, tpi^dh, where dh is hole creation operator. 
The first order correction of expansion has the following form 



S-+ cx < ijR{x, t)ijR{x - e, t)dl{x - e, t) cos[kFX + V^4>{x, t)]e^['^«'^«("'*^+^^'^"^("'*)l 



^^hnVu{o.o)+^L'PLio,o)] cog[y^^(o^ o)]4(e, 0)Vilj(e, 0)^Jj(0, 0) > e'(2'^---^"+^)- 
Similar to Eq. p.8p we could decompose S ^ to two terms as follow 

S-+{x,t) = e'^''-'"^+'"'^^I^{x,t) + e'^''-'"'+^'"'^''I+{x,t) (4.19) 
where the expression for is similar to Eq. (|3.9p with exponents given by 

2 

, (0 -V- ^\^^K^ 

4 \^:^/K ^/K 



R 

y± - 

^ 4 V7^^/i^ ' ^/K 



Vf = -\^ + ^- {2±l)^K\ (4.20) 



In general, the higher order correction will include more powers of V'ljV'i a-^d general expression would be 

4")± = if^ + ^ + (2n±l)Vif'' 



-{n)± If IL 



2 \^ 
--j=-{2n±\)'jK\ (4.21) 



20 



with momentum given by 

A:^") = TT - kh + i2n+l)kF 
and the energy of the excitation for given momentum q is 

w„(g) ±e[-q + {2n + l)kF + tt]. 



(4.22) 



(4.23) 



ujn and u!„ for n = 1, are depicted by dashed hues in Fig®, which represents vanishing singularities of spectral 
function. The sum of exponents is 



, -(n) 1 
4 + = 2 



(271+ 1)^^^ 



IR - IL 

2ttJ~K 



}_ flL+ 7-R 

2 



(4.24) 



We see that our results Eq. (|4.1ip and (|4.24p are the same as Eq. (17) of Ref [1], which is the singularity exponent 
for boson creation operator, upon identifying 8± 1r/l/^/K. The actual values of these phase shift parameters are 
in general different in the two models however, being determined by Galilean invariance in the Bose gas model. This 
correspondence might have been anticipated since a boson creation operator is related to the corresponding fermion 
one by a Jordan- Wigner string operator— just as is the S*^ operator. Furthermore, the xxz model is equivalent to a 
lattice boson model with an infinite on-site repulsion which restricts the occupancy to or 1. 

Now we look at the effect of Umklapp scattering term at zero magnetic field. The Umklapp scattering term is in 
the following form 



'^Ri^)'^R{x)i'L{x)ipL{x) +h.c 



(4.25) 



At non-zero magnetic field we could ignore this term as it is highly oscillatory, due to e*""°^^ prefactor. At zero field 
we have Akp = 2tt; thus this term does not oscillate and we need a more careful treatment. In the bosonized form of 
the XXZ model, it can be seen that the Umklapp term is irrelevant for < A < 1, and is marginal for A = 1, thus 
in this regime we could look at the effect of this term perturbatively. We focus on the effect of this term on S ^. In 
general the higher order Umklapp term could be written as follow 



m-t-m 



2m+m 



(4.26) 



Where m and ffi are arbitrary integers, and this term is actually 2m -|- m order in the Umklapp perturbation. These 
operators gives us zero corrections unless we keep higher powers of the m{x) term, Eq. |4.1| . at least to power 2m. 
Thus the most general non-zero term of both Umklapp and oscillatory term to(x) is in the following form 

Pn^m,rn = C/am+m (V^l (.t) Vi?, (^) )™ (^I (0) V-i?. (0))" (4.27) 
Now by bosonizing the above expression and plugging it into the definition of S ^ we have 

S-+0, J < dix) cos(fcj. + V^(^(.T))e7SK(^'«'^«(-)+^'^'^^(-))e»^™(™-")v(-)e-*2^/4;^7^E™+* vi^.) 

ill rn+m 

X e'2v^E'^^to0cos(V^(^(0))e^^^«'^«^°)+^"'^"(°»e^^^("+")^(°)dt(0)>nrf'y« "Q ^2^^ ^^ ^S) 

The effect of this term on the singularity exponents of S ^ could be evaluated by power counting and the result is 
given by the following expression 
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+ K{m^ + 2m + Afh) - (m + 2m) 



(4.29) 
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The — (m + 2m) term in above equations conies from the integration variables Zj,yi. The overall exponent is given 
by the summation of these two exponents, thus we have 

As we are considering the effect of Umklapp term at zero field, Akp = 27r; these are singularity exponents of the 
spectral function at threshold frequencies given by Eq. |4.9| . Compared to Eq. [4.11| . we see that for fixed n 
higher order Umklapp terms give larger and larger exponents, and are more irrelevant. Therefore, at each threshold 
energy, we get a set of singularities with progressively weaker exponents. In this case, these corrections seem truly 
unimportant. That is, we don't get any new singular energies, just sub-dominant corrections to the singularities at 
the energies we already have. 

One important point is that if we take n,m — but arbitrary m, at half filling and for A = 1 wc have 2K — 1 = 0; 
thus higher order non-chiral Umklapp corrections do not change the singularity exponents, based on power counting. 
But wc should be careful at that limit, because our result was based on power counting; in general upon evaluating 
the integrals more carefully, there could be some logarithmic corrections to the correlations functions which could 
change the behaviour of spectral functions near singular frequencies. 

Till now wc only considered the effect of Umklapp terms only at zero magnetic field. At finite magnetic field the 
Umklapp interaction has the following form 

Uu ^ gu cos4{V^(p{x) + (kp - TT/2a)x) (4.31) 

Where a is the lattice spacing. In general this term is oscillatory and could be drop out at low energies. But for 
weak enough magnetic field, the wavelength of the oscillation is very long, therefore in that limit this term should 
be treated carefully. We claim that at low enough temperatures and weak magnetic field the Umklapp term affects 
neither the threshold frequency nor the singularity exponents, but it will change the overall behaviour of the spectral 
function; the reason is as follow. If we include the effect of Umklapp term perturbativley, it could be easily shown 
that such higher order terms can not change the oscillation wave-vector of the spectral functions. It only modifies 
the non-oscillatory part of the spectral function without changing the oscillatory part. Therefore, if the oscillations 
wave-vector remains intact the threshold frequency does so. 

The singularity exponents does not change because, to find the singularity exponents of spectral functions, we need 
to study the behaviour of spectral functions at frequencies, w, around the threshold frequencies, w^. In principle the 
probe frequency could be chosen as close as possible to the singular frequency such that the Umklapp term effects 
would be irrelevant at those energy difference scales, \u — uii^\ <^ {{kp — 27r/a)w|. Thus the Umklapp term would not 
change the singular exponent at ujl', we expect a cross over regime where the effect of Umklapp will be important at 
energies near to |a; — lul \ sa {kp — 27r/a)v\. 

V. ELECTRON SPIN RESONANCE WITH UNIFORM DZYALOSHINSKII-MORIYA INTERACTIONS 

Electron spin resonance provides a sensitive probe of spin dynamics. A microwave field is weakly Zeeman coupled to 
the (7 = components of the spin operators. In the standard (Faraday) configuration, the microwave field is polarized 
perpendicular to a static magnetic field. For simplicity we restrict ourselves to the relatively simple situation of 
Eq. (|l.ip . with DM vector and magnetic field in the z-direction. Then, as discussed in Sec. I, a uniform DM 
interaction added to the xxz model of Eq. (jl.ip simply shifts the parameters J and A and the momentum, q in 
the transverse spectral function. Therefore the ESR adsorption intensity is proportional to the transverse spectral 
function at q = a = arctan(£)/ J). Low temperature ESR measurements on quasi-lD antiferromagnets with uniform 
DM interactions could therefore probe the edge singularities predicted by X-ray edge methods that disagree with 
standard bosonization results due to the effects of band curvature. Such ESR results would be especially useful if 
they were done with circularly polarized microwave radiation since then S and could be measured separately. 
A further major challenge for such ESR experiments would be that the theoretical predictions give S'^^ {a,u!,h) for 
fixed /i as a function of tu. However, in an ESR experiment, uj is normally fixed at the resonant frequency of a 
microwave cavity and h is varied, w can only be varied by using a sequence of microwave cavities with different 
resonant frequencies. Alternatively, theoretical line-shapes could be produced for fixed w and varying h but these 
would be complicated since the critical exponent a varies with h. For simplicity, wc just discuss the line shape versus 
frequency at fixed h here. 
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We begin by discussing the T = limit. The ESR adsorption intensity can be simply read off from the results of 
Sec. II, III. We first consider the case of zero static field with the microwave field in the xy plane. Then the adsorption 
intensity has a lower threshold near which: 

9{uj - wsina) 

« 7 ■■ xTT (^-1) 

[uj — V sina)^' 

with n = 2 — 1/{2K) — K. Here v and K are determined in terms of Ae// = A cos(q;) by Eq. (|2.10p . For small a we 
expect the results of standard bosonization to apply at somewhat higher energies, w — fa 3> a'^ J. In this region we 
obtain: 

/ a; cx ) '-—^ 5.2 

with 77 = (1 — 2ii')^/ (8/-r). Note that 1 — 2?^ = 2 — \/{AK) — K ^ ji, a. different exponent than occurs at the threshold. 
As discussed in sub-section IIA, /(cj) in Eq. (j5.2p is non-monotonic, eventually passing through a minimum and 
starting to increase again as oj increases. However, since the formula is only valid in the low energy regime, uj J, 
whether or not this minimum occurs in the frequency region where the formula is valid depends on a and Ae//. If 
Ae// > cos a, the value resulting from A = 1, then the minimum predicted by standard bosonization is not in the 
region where the approximation is valid. In this case the intensity is monotone decreasing up to high frequencies 
where our techniques break down. Ae// is typically close to 1. In fact, with some assumptions about the higher 
energy levels of the magnetic ion, it is exactly one<^ In this case, the transverse spectral function becomes the same 
as the longitudinal one discussed extensively in @, [l^ Hi], and[33|. The edge exponent has the value /i = 1/2, first 
obtained from the 2-spinon approximatioi*^ in this case and the standard bosonization prediction of Eq. (j5.2p . with 
r] = 0, reduces to a (5-function which fails to capture many features of the actual spectral function for non-zero a. In 
particular there is a narrow peak of width cx followed by a slowly decaying tail at higher uj. 

At non-zero field, h, we may again use the results of Sec. II, HI, which are less complete in this case. The X-ray 
edge results of Sec. HI imply a quantum phase transition as the magnetic field is increased, occurring when the 
field-induced magnetization, m{h), obeys a ^ H = 2T:m w Kh/{Ti/v) as can be seen from Fig. [6]. S has a 
threshold singularity at a frequency lol ~ v{a — H) for H < a <^ 1. A threshold singularity was also predicted in 
HI for S'^ in this field range, at a higher frequency, of approximately v[a + H). However, as discussed at the end 
of that sub-section, we expect this to be broadened. On the other hand, for H > a, we expect to have a sharp 
2-sided singularity at a frequency of approximately v{H — a). The other threshold singularity predicted for at 
v{H + a) in IIB is likely to be broadened. (The presence of two peaks at these energies was first predicted in [s^ 
and was observed experimentally in [20|.) The precise energies and critical exponents for these singularities could be 
predicted by numerical Bethe ansatz calculations but analytic expressions are not available. Assuming a, iJ <C 1, 
we expect the spectral functions to cross over to the form predicted by standard bosonization at energies somewhat 
higher than the singularities, {H — a^aJ <^ v\H — a\ <^ J: 

As discussed above and sketched in Fig. [1] , these functions are non-monotonic but the minimum only occurs in the 
energy region w ^ J, where the approximation holds, for a certain parameter range of a, A and h. In addition to 
these dominant singularities, as discussed in IIIB, we expect many additional weaker vanishing singularities extending 
down to low energies. 



Finite Temperature Broadening 

At finite temperature the sharp peaks are broadened. This can be calculated using standard bosonization for 
{H - aYaJ <t: J. At finite T the spectral function of Eq. becomes: 

/ cx Im [sin(27r77)(27rr)''B(l + ry - i{uj + va)/{ATTT),-l - 2i^)B{i] - i{uj - va)/{ATTT), 1 - 2??)] . (5.4) 

where B(x, y) = T{x)T{y) /T{x + y) is the Euler beta function. For weak anisotropy, small a and 1 — ^eff, and hence 
?/ <C 1, a Lorentzian line-shape occurs: 



[uj — va)"^ + (47rr?/)2 



23 



with a similar broadening at finite H. The width of the peak is AttTt] « 2T(1 — Ae//)/7r for Ae// close to 1. This 
is essentially the same result derived in [l^ for the ESR width due to exchange anisotropy parallel to the magnetic 
field. 

In section imi we found the singularity exponent and behaviour of the transverse spectral function near the thresholds 
at zero temperature using X-ray edge methods. Now we look at the finite temperature effects on the spectral functions, 
and we find that finite temperature results in a non-Lorentzian broadening of the transverse spectral functions near 
the thresholds for \uj — e{q)\, T ^ /2m. Non-zero temperature has two effects on the spectral function Eq. (|3.2ip . 
First, the Green's function for excitations near the fermi surface, S^^'^ is modified to the usual finite T form by a 
conformal transformation. Secondly, the Green's function for the d operators would be modified to finite T form 
- with step functions 6{e) replaced by Fermi functions nj7(e). This would allow a hole contribution to S and a 
particle contribution to S'^ . However, at low T <^ t{kp) these would be negligible. So the only important effect may 
be in 5"^°' . Thus in this regime we have 

S-+(q,u})(x dxdt ^— ^-^ '- -, (2kF - tt < q < n) 5.6 

J {sm{2nT{e + i{t- x/v)))Yi{sm{2TTT{e + i{t + x/v)))yi 

By doing the integration over x we have 

S~+(q,uj)(x / dt '- i ^ ^, {2kF - 7T < q < tt) (5.7) 

J-oo (sin(27rT(e + i{l - u/v)t))y4 (sin(27rT(e + i{l + u/v)t)))''L 

By evaluating the above integral we can get the finite T behaviour of the spectral function. From the above equation 
it can readily be seen that the spectral function is a pure real number; in the appendix we will prove that it is positive 
too. Let us look at the behaviour of above spectral function for 2kp — tt < q < tt, at h = and Ag// ~ 1, small 
anisotropy 77 1 so that vr = 1/2 — 77 and = 0. In this regime we have 



^ Re 



(5.8) 



In Fig. [7] we have depicted, the spectral function for different values of temperature. As temperature gets higher the 
broadening increases; as can be seen, for small T the broadening is asymmetric and so is non-Lorcntzian. Furthermore, 
the width is 0{T), not suppressed by a factor of 1 — ^eff as predicted by standard bosonization. 

Since the effects predicted by the new theoretical methods occur at very low energy scales, a highly one-dimensional 
spin compound would probably be needed to observe them, in order that three dimensional exchange processes would 
be negligible. Furthermore, materials like KCuGaFg with both uniform and staggered intra-chain DM interactions 
may not be suitable since the staggered DM interactions tend to have a larger effect than the uniform onesJ^ Thus 
identifying the right material to test these predictions remains an open chall eng e. An alternative approach might be 
to study metallic quantum wires with spin-orbit couplings. It was shown in [37| that analogous phenomena occur in 
that system. 



VI. CONCLUSIONS 



By applying X-ray edge techniques, we have obtained results on the transverse spectral function of the xxz spin 
chain in a magnetic field. We illustrated why standard bosonization techniques fail near threshold energies, even 
when these occur at low energy. In the zero field case we have exactly determined a critical exponent governing the 
lower edge singularity, for all |A| < 1 and all wave- vector. For the finite field case, we have shown how this exponent 
can be determined from parameters which can be obtained from solving Bethe ansatz equations and which also 
determine the behaviour of the longitudinal structure function, fermion spectral function and the finite size spectrum. 
We have argued that, for general magnetization, a large number of increasingly weaker singularities occur in the 
spectral function, extending all the way down to zero energy. We derived results for the finite temperature spectral 
function using X-ray edge methods, obtaining strikingly different behaviour than that given by standard bosonization 
at < 1 — Ae// ^ 1. The line-shape is non-Lorcntzian and the line width is 0(r), unsuppressed by 1 — ^eff- We 
pointed out that electron spin resonance measurements on spin chain compounds with uniform Dzyaloshinskii-Moriya 
interactions would provide a way of experimentally confirming, for the first time, the new bosonization results being 
obtained on spin chains, using X-ray edge techniques. 

In the case of staggered DM interactions, the most interesting ESR signal occurs when the magnetic field is transverse 
to the DM vector. This may also be the case for uniform DM interactions, but we leave this for future work. 

We thank J-S Caux, Ion Garate, Leonid Glazman, Rodrigo Pcrcira, Olcg Starykh, Hidekazu Tanaka and Izumi 
Umegaki for helpful discussions. This research was supported by NSERC and ClfAR. 
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FIG. 7: The behaviour of 5+ (w,q) speetral function for at finite temperature for the case i^r ~ 1/2, i^l = 0, 
corresponding to A = 1, /i = 0. We see that for smaU enough T the broadening is asymmetric. 



Appendix A: Positivity of spectral function 

In this appendix we will prove that spectral function Eq. (j5.7p is positive. We have 

g4[w-e(fcp)]i^27rT)'^«+'^i' 



S oc 



dt 
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CJO 

oo 



POO /• 

L 



(sin(27rr(e + i{\ - u/w)t)))''«(sin(27rr(e + i{l + u/v)t))YL 
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(sin(27rT(e + i(l - u/v)t))Y'i{&m{2'KT{e + i(l + u/v)t))Y^ 



(Al) 



By defining the branch cut on (— icx3, — ie) lj(ie, ioo) and doing the change of variable i — > — i in second integral, we 
have 



S '^{q,Lj) oc Re 



exp —i—yvr + sgn(l — u/invn) / dt- — , , ^ , — ; — , , — 

2^ ' J Jo I sinh((l-M/w)t)|''«|sinh((l + 



(A2) 



We will prove shortly that above integral is positive in general, but let us first look at the special cases of either 
vr = or = 0^ which is relevant for weak anisotropy at zero magnetic field, so m < t;. 
Suppose that 1/^ = thus we have 

gi[w-£(fcp)]t/(27i-T) 



S'-+((j,w) OC (27rT) 



Re 



dt- 



oc {2nT)'"'~^Re f exp (-i^i^i?) B 
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(A3) 



Where B\x, lA is Euler's beta function. Let's define W _ . 

L ' -^-1 47rT(l — u/u) 

terms of gamma function we have 



47rr(l - ujv) 

— — and by using the definition of beta function in 



S-+{q, uj) oc {2nTy-'Re exp -i-i^R 



T[-iW + iyR/2]r[l - i^r] 



(A4) 



2'- J T[l - {iW + ur/2)] 
Then by using the identity that r[l — z]r[z] = 7r/sin(7rz) we have 

S-+{q,oj) oc {2nTY^-^r[~iW + VRl2\V{iW + J^ii/2]r[l - VR\Re (exp sin7r(iVF + vrI2')^ (A5) 
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by "writing the sine function in exponential form, we finally get 



(A6) 



so we see the spectral weight is positive, for all values of vr. This expression also holds for the case oi u > v. 
Actually what we have shown is that the Fourier transform of 1/ sin(27rT(e + i{l — u/v)t)Y is given by a real positive 
function, in the following form 
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Where F.T stands for Fourier Transform. Now by using above equation and taking the convolution of Eq. (|5.8p . we 
have 
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We see that all the functions in above equation are positive; thus whole the integral is positive. 
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